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Abstract 

We propose a new method based on the invariant properties of the ansatz of 
the Hirota method whcih have been discovered recently. This method allows one 
to construct new solutions for a certain class the dissipative equations classified 
by degrees of homogeneity. This algorithm is similar to the method of "dressing" 
the solutions of integrable equations. A class of new solutions is constructed. It 
is proved that all known exact solutions of the FitzHygh-Nagumo-Semenov equa- 
tion can be expressed in terms of solutions of the linear parabolic equation. This 
method is compared with the Miura transforms in the theory of Kortveg de Vris 
equations. This method allows on to create a package by using the methods of 
computer algebra. 
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We consider the quasilinear parabolic equation 

(1 + biu + b2V^)ut - {hi + h2U + h:iu'^)ux - goiu^Y 

4 

-{kQ + kiu)uxx + ^'^iU^ = ^- (1-1) 

i=l 

where bj{x,t), hj{x,t), kj^i{x,t), j = 1,2, ipi{x,t), i = 1,...,4, gQ{x,t) are smooth 
functions or constants. 

It is assumed that the coefficients of Ut anduxx do not vanish in the range of the 
function u. 

We seek a solution in the form of the fraction (R. Hirota): 

(1.2) 

After multiplication by F{x,t)'^, for the functions G{x,t) and F{x,t) we receive a homo- 
geneous equation of the fourth order of homogeneity with respect to all functions. 
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A special case of Eq. (IJ.) has the form 



(1 + biu)ut - {hi + h2u)ur, - fcoMxx + V^i"* = 0' (1-3) 

i=l 

which, after the substitution of ( p. .21) and multiphcation by F(x,t)'^ becomes an equation 
of the third order of homogeneity with respect to the functions G{x, t) and F{x, t). 

Thus this equation is of the third order of homogeneity, and the the components or 
(ux)"^ comphcate the equation so that it is an equation of the fourth order of homogeneity. 

Special cases of this equation are widely known (see the extensive bibliography in 
I, |,0,[Tl,[ll). 

If hi = 0, (fi = 0, i = 1, 3, then this is the Biirgers equation in the wave theory [Q. 

If hi = 0, /i2 = 0, ipi = —1, !f2 = —fi, f3 = 0, then this is the Fisher-Kolmogorov- 
Petrovskii-Piskynov (FKPP) equation (see the bibliography in 0, |^, |13|); 

If hi = 0, /i2 = 0, ipi = —1, ip2 = 0, (/93 = —ipi, then this is the FitzHygh-Nagumo- 
Semenov (FhNS) equation in the wave theory (see bibliography in HTj-^]), fT^i P- 12), 
and the Allen-Cahn equation (which is close to the FhNS equation) in the theory of 
phase transitions |]TD|] and in ecology [0, |T3]. Equations in which the coefficients of the 
Laplace operator depend on the function u appear in models of the theory of epidemic 



distribution |20|. Equations in which the coefficients of the derivative Ut depending on u-^ 



arise, for example, in models of processes of ocean water freshening |2T|]. The initial 



reasons on a theme of the given operation can be detected in operation |2^ . 

The method proposed for constructing solutions of Eq. (1.1) ideologically goes back 
to the Hirota method. An essential distinguishing feature of the present paper is that a 
new invariance property of solutions of equations from a given class is discovered (cf. [16, 
17]). The equations obtained from (1.1) after the substitution of (1.2) can be arranged in 
the order of homogeneity as embedded equations. 

From now on, we use the following notation for the derivatives: 

d,{-) = {■),, u = x,t (1.4) 

We choose the solution (1.2) of Eq. (1.3) in the form 

G{x,t) = exp{H{x,t))B{x,t) + Q{x,t), (1.5) 
F{x,t) = (l + A{x,t)exp{H{x,t))y (1.6) 

All the other representations of the solution obtained by the Hirota method can be 
reduced to this form. There is only one possibility to extend this notation, namely, to add 
sums of series in powers of the exponentials exp{iH{x,t)), i = 2, . . .. Let us substitute 
(|1.5|) first into ([L^) and then into ( |1.3|) and equate the coefficients of equal powers of 
exp{H{x, t)) (these powers are 0, 1, 2, 3) with zero. 

We denote the obtained system of equations by (egO), (egl), (eg2), (eg3). 

Equation (egO) is the same original equation ( p..3|) for the function Q{x,t). Hence 
Q{x,t) is a solution of this equation. 

Let us consider the equation obtained at a power of exp{3H{x,t)), namely, Eq. (eg3): 

^(fiVs - 2koiA,f + B{-biAt + 

+A[b\^i + y?3) + 2koA^B^ + B{-At + BiBt + hA^ - h2B^ + k^A^ 
+A\Bvi + Bt- hiB., - KoB^^) = 0. (1.7) 
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The special case considered in [|], p. 51, Section 24, and in [2], p. 190, is the case of 
B = (±)v4, where the equation has the form 

- A^B + B^ - ABAt - 2BAl + AB A^^ + A^Bt + 2y4A^.5^ - A^B^^ = 0. (1.8) 

and is an identity. 

First, we present a lemma in which we prove that, in the conventional ansatz form 
of a solution, there is considerable arbitrariness related to the "extra" function A{x,t), 
which must be excluded. Namely, this function must be added to the function H{x, t) 
according to the formula H{x, ^)(new) = H{x, t)j^Qj^^+exp(ln(A)). Then the same notation 
^))(new) is preserved for the function H{x,t). 

Lemma 1.1 Suppose that the functions B{x,t), A{x,t) in (1.5) belong to C'^[R^ ® R^] 
and satisfy Eq. (eqS), and the function B{x,t) has the form 

B{x,t) = M{x,t)A{x,t). (1.9) 

Then the function M{x,t) satisfies Eq. (1.3) 

3 

(l + h{x, t)M)Mt - (/ii(x, t) + /i2(x, t)M)M^ - ko{x, t)M,, + ^ ^i{x, t)M' = (1.10) 

i=l 

for all functions A{x, t). 

Proof. The proof of the lemma is directly obtained by the substitution of (|1.9|) into 
( |1 . 7|) . Equation (1.10), which follows from Eq. (1.7), coincides with (1.3). Hence the 
function M{x, t) is a solution of this equation. Below M{x, t) stands for some known 
exact solution of Eq. (1.3). □ 

It follows from Lemma 1.1 that the number of arbitrary functions in the ansatz can 
be decreased. Namely, we write A = exp(ln(A)) and add this summand to the function 
H{x,t). In the present paper, solutions of the form ( |1.2|) are called solutions of the first, 
second, etc. mode with respect to their complexity (the number of terms containing a 
power of exp(iif) as a factor). 

Definition 1.1 The set Q of solutions of Eq. (1.3) of the form 

exp{H{x,t))M{x,t)+Q{x,t) 
U[X,t) = -— , (l-ll) 

where M {x , t) , Q{x,t) are some solutions of Eqs. (|1.3| ) and H{x,t) is an arbitrary func- 
tion, will be called SOLUTIONS OF THE FIRST MODE; solutions of the form 

^ expjHjx, t))M{x, t) + exp(2g(x, t))R{x, t) + Q{x, t) 
' l + Zi(x,t)exp(i/(x,t))+exp(2/f(x,t))Z2(x,t) ' ^ ' ' 

will be called SOLUTIONS OF THE SECOND MODE; etc. 
We have the following theorem. 
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Theorem 1.1 Suppose that all coefficients in Eq. (|1.3| ) are constants. The function M 
is a solution of the equation 



(1 + h^M)Mt - ih + h2M)M, - koM,, + ^ = 0, (1.13) 

1=1 

and Q is a solution of the equation 

3 

(1 + hQ)Qt - (/ii + /i2Q)Q. - fcoQx. + E v>Q' = 0, (1.14) 

i=l 

In addition, they are related by the compatibility relations 
{M - Q)b^koH^^{x,t) 

+H^{x, t)(M{bihi - h^) + Q{-bihi + h^) + 2biko{M^ - Q, 



+ {2 + Mb^+Qbi)koH,{x,t)^ 

+{M - Q) (^M%ip3 - Q%iip3 + M{-ip3 + bi{ipi + ip3)) - Q{-ip3 + bi{ipi + (^3)) 

+b,{b,{Mt - Qt) + h2{-M^ + g,.)))) = 0, (1.15) 

and, together with the function H{x, t) satisfy the equation 

- (M - Q)b^Ht{x, t) + H,{x, t)(M/i2 - Qh2 - 2koH,{x, t)) + (M - Q) Vs = 0. (1.16) 
Then the function u{x,t) is a solution of Eq. ( |1.3| ). 

(The arguments of the functions are omitted here for brevity.). 



Remark 1.1 The equation ( |1.13|) is the Riccaty equation for the functions H^-^'^\x,t), 



and equation ( |1.14|) is the Hamihon-Jacobi equation for the function H{x,t). 



If bi = 0, then the statement of Theorem 1.1 becomes simpler and the function H{x, t) 
can explicitly be written in terms of the functions M{x,t) and Q{x,t). 



Theorem 1.2 Let bi = 0. Suppose that all coefficients in Eq. ( |1.3|) are constants. The 
function M is a solution of the equation 

3 

Mt - {h, + h2M)M, - fcoM,, + E V^M' = 0, (1.17) 

i=l 

and Q is a solution of the equation 

3 

Qt - {h, + h2Q)Q. - koQ,, + E y^iQ' = 0. (1-18) 

i=l 

In addition, they are related by the compatibility relations 

{Mt - Qt) dx^ih - m) - 2(M2 - Q^){hl + Uko^s - /12m) 
-4(M - Q){hih2 + Ako^2 - him) + 16A;oC^(i) 

+4fco(/i2 + 3m)(M,-Q,) = 0, (1.19) 



where m = ((±)y^ (h"^ + Skoips)), and the function H{x,t) has the form 



Ako 



Then the function u{x^t) is a solution of Eq. (|1.3|) . 

Proof. We give a brief proof of this theorem. We substitute ( p..ll| ) into ( |1.3| ), reduce 
to the common denominator, and consider the numerator of the obtained large equation. 

We equate with zero the coefficient of the exponential raised to the zero power, i.e., 
the coefficient of exp(O). The equation thus obtained is exactly Eq. (|1.3| ) for the function 
Q{x,t). Hence the function Q{x,t) is a solution of this equation. 

We equate with zero the coefficient of the exponential exp(3if(x, t)). The equation 
thus obtained will be called Eq. (eg3). Again it is exactly Eq. (|1.3|) . Hence the function 
M{x, t) is a solution of this equation. 

We equate with zero the coefficient of the exponential exp{H{x, t)). The equation thus 
obtained will be called Eq. (egl). 

We equate with zero the coefficient of the exponential exp{2H{x,t)). The equation 
thus obtained will be called Eq. (eg2). Since the functions Q and M satisfy Eq. (1.3), we 
exclude the second derivatives of the functions Q^^ and M^^ from Eqs. (egl) and (eg2). 

We express H^x from Eq. (egl) and stubstitute it into Eq. (eg2). Then we obtain 
the Hamilton-Jacobi equation (1.16). For bi = this equation in Theorem 1.2 naturally 
determines the function H{x,t) (1.20). If we express the derivative Ht from Eq. (1.16) 
and substitute it into Eq. (egl), then we obtain the compatibility condition (1.15) or 
(1.19), respectively. If we substitute the derivative Ht into Eq. (eg2), then we also obtain 
the compatibility condition (i.e., the two equations coincide). In conclusion we note that 
there exist two branches in the relations given in the theorem with different signs. The 
proof of the theorem is complete. □ 

In fact, the theorems proved above describe the procedure of constructing new solu- 
tions of the equations (an analog of the superposition of solutions), which is specified by 
the properties of Eq. (|T]3|) and the ansatz. 

Let us consider two smooth solutions of Eq. ( |1.3|) : M{x,t) and Q{x,t). The following 
function is assign to these solutions: 

_ Q{x,t) + M{x,t)exp{H{x,t)) 

l + exp{H{x,t)) ' ^ 

where H{x,t) is determined by (1.16) or by formula (1.20) if bi = 0. 
The following theorem states that the solutions are invariant. 

Theorem 1.3 Suppose that two pairs of functions 

(Mo{x,t),Qo{x,t)) and (^Mi{x,t),Qi{x,t)) 

are solutions of the equations 

3 

(1 + biMj)Mj^ - {hi + h2Mj)M,^ - koM,^^ + ^iM/ = 0, (1.22) 

i=l 
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3 

(1 + biQj)Qj^ - {hi + h2Qj)Qj^ - koQj^^ + 'PiQ/ = 0, (1-23) 

i=l 

where j = 0, 1, and that they determine the solution by formula (1.11). Suppose that the 
assumptions of Theorems 1.1 and 1.2 hold for each pair. Then the two functions 

exp{Hi{x,t))Mi{x,t) + Qi{x,t) 



Q{x,t) 



l + exp{Hi{x,t)) 



^ exp{Ho{x, t))Mo{x, t) + Qo{x, t) 
1 + exp{Ho{x,t)) 

are solutions of the equations 

3 

(1 + hiM)Mt - {hi + h2M)M^ - koM^^ + ^ ^iM' = 0, (1.25) 

i=l 
3 

(1 + biQ)Qt - {hi + h2Q)Q^ - koQ,, + E ¥^iQ' = 0, (1-26) 

i=l 

and the function 

_ exp{H{x,t))M{x,t) + Q{x,t) 

l + exp{H{x,t)) ^ 

is a solution of Eq. ( |1.3|) . 

Proof. Each of the functions Mj, Qj, j = 0,1,, satisfies Eq. ( |1.3|) written for this 
function. Pariwise, they satisfy the corresponding compatibihty condition ( |1.15| ) or ( |1.19| ). 
The theorem can be proved by a direct substitution. □ 



Remark 1.2 After calculating the derivative with respect to x, it is possible to consider 
the compatibility condition (|1.15|) as a modified Biirgers equation with a potential. 

Corollary 1.1 The class of solutions (|1.271) is closed with respect to the procedure de- 
scribed in Theorem 1.3. Thus ( |1.27|) defines a surjective mapping: 

Vol: n. (1.28) 

Similar results take place for Eq. (1.1). We perform the change of the variables = 
0201, ^2 = —(01+02+0102), V^3 = l+oi+a2, = —1, which is equivalent to a source-sink 
function written so that the four roots are explicitly distinguished: u{ai—u){a2 — u){l — u). 

We have a lemma similar to Lemma 1.1 and the following theorem. 

Theorem 1.4 Suppose that the solution (1.1) has the form (1.11) and the relations hold 
in one of the following cases: 

1) {Q{x,t) or M{x,t)) = const = —ko/ki, {ko + ki){ko + aiki){ko + 2aiki) = 0; 

2) 61 = ki/ko, 62 = 0; 

3) bi = 0, 62 = 0, ki = 0. 

Then the functions M andQ are also solutions of an equation coinciding with Eq. (1.1). 
Moreover, these functions are related by the compatibility relations {the Riccati equation 
for Hx{x,t)) . In addition, the Hamilton- Jacobi equation must be satisfied. 

(Here we do not write these equations, since they are very cumbersome.) The proof 
is similar to the proofs of the preceding theorems. 
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In the case of variable coefficients, there are additional possibilities to satisfy the system 
of equations. Since we construct exact solutions, not all coefficients of the equation are 
arbitrary. 

We have the following theorem. 

Theorem 2.1 Suppose that bi{x, t) = 0, the function M{x, t) is a solution of the equation 

3 

Mt - (hi{x,t) + h2{x,t)M)M, - ko{x,t)M,, + Y,^i{x,t)M' = 0, (2.29) 

i=l 

the function Q{x,t) satisfies the modified equation 

6 

Qt - {hoi{x,t) + ho2{x,t)Q)Q, - koo{x,t)Q,, + J2^i{x,t)Q'-'' = 0, (2.30) 

and the coefficients in the equations satisfy the relations 

Q(^(pi{x, t) - <^^{x, t) + Q{ip2{x, t) - ip^{x, t)) + (5^((^3(a;, t) - ipQ{x, t))) 
+ {hQi{x,t) - hi{x,t) + {ho2{x,t) - h2{x,t))QjQx 

-{koQxx - KoQxx) = 0, (2.31) 
and the functions M , Q and the coefficients satisfy the second compatibility equation 

+2(M - Q){±hl{M^ - g2) + 2hi{M - Q){{±)m'^ + /lam) 

+/i2'((M - Q)(M + Q)m - i±)2Mko, + i±)2Qko^ + (±)6A;o(M, - Q,)) 

+2h2{{±)4ko^3{M' - Q') + {±)koh2,{M -Q)- m{k^,{M - Q) + ko{M, - Qx)) 

+2ko(m(2{M - Q){2^2 + ^Vz{M + Q)) - 4C^(t) + /i2,(M - Q) 

-(±)4(^3A;o.(M - g) + (±)24v93A;o(M, - Q,) + (±)4^3,A;o(M ' Q))) = 0, (2.32) 

where m = {\J (/i^ + S/^oV^s)) and the function H{x, t) has the form 



r (M(x, t) - Qix, t))( - h2{x, t) ± Jh2{x, tf + 8A;o(x, t)^^{x, t)) 
Hix, t) = Coit) - J dx. 

(2.33) 

Then the function u{x,t) (|1.11|) is a solution of Eq. (|1.3|) . 

Proof. The proof is similar to that of Theorem 1.2. Note that, ffist, in this case the 
set of solutions to the system of equations for the functions M, Q and the coefficients of 
the equation described in the theorem is not empty; 

second, it is possible to construct solutions with variable roots of the source-sink 
function. 
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We substitute (1.11) into (1.3), reduce to the common denominator, and consider 
the numerator of the obtained equation. By equating with zero the coefficient of the 
exponential raised to the zero power, we obtain Eq. (egO), which is the compatibility 
condition for coefficients (2.22). Prom this condition we express some coefficient and 
exclude this coefficient from the other equations. The equation obtained by equating the 
coefficient of the exponential cxp{H{x, t)) with zero is called (egl). The equation obtained 
by equating the coefficient of the exponential exp{2H{x, t)) with zero is called (eg2). Since 
the function Q satisfies Eq. (2.30), we exclude the derivative Qt{x,t) from Eqs. (egl) and 
(eg2). We can express Hxx{x,t) from Eq. (eql) and substitute it into Eq. (eg2). Then 
we obtain the Hamilton-Jacobi equation, which naturally determines the function H(x, t) 
(2.24). The proof of the theorem is complete. 

In the construction of the solution of the second mode, the transformation relates the 
solution of the original equation and the solutions of two modified equations, and so on. 

The FitzHygh-Nagumo-Semenov equation 

EXAMPLE 1. We apply our method to the FitzHygh-Nagumo-Semenov (FhNS) 
equation known in wave theory: 

ut - ^iU^x/ (2a^) - ip2.u + ip'iU^. (2.34) 

where ko = fs/ (2a^) and (pi = (p^ are constants. We show that this method leads to new 
results and provides a new understanding of the results (see [3], p. 64, and [9]). 

First we show the action of the semigroup related to the transformation (1.11) and 
explain how the translation constants pass into new solutions, by choosing the form of the 
functions (of its solutions, [1], p. 17, [3], p. 51) satisfying the assumptions of Theorem 1.2. 
Namely, we consider the functions u{x,t) — M{x,t) or u{x,t) = Q{x,t): 

M{x,t) = -{1 + ex.p{mo + ax + bt))~^ , (2.35) 
Q{x,t) = (-l + exp(2ax + go))/(l + exp(2ax + go)), 6 = -3(/?3/2. 



We calculate the function 



^ l + exp(mo + aa;-3W2) 
^ ' 1 + exp(go + 2ax) ^ ' 

and obtain the solution 

^ l-exp(-go-2ax + ln(2)) 

1 + exp(— go — 2ax + ln(2)) + exp(mo — go — ax — 3t(p3/2) 

Thus the original transition constants go, uiq pass into a new solution and a new constant 
In (2) appears. 

Note that this is a well-known solution that describes the wave interaction ([2], p. 190, 
[3], p. 51, and [7]). Below we prove that this solution can be expressed via the solution 
of the linear parabolic equation. 

There exists another version of the functions for which the assumptions of Theorem 1.2 
are satisfied: 

M{x,t) = {l + exp{mi- ax + bt))-^, (2.38) 
Q{x,t) — (1 — exp(gi — 2ax))/(l + exp(gi — 2aa;)), b — —3if2/2. 
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Then the function H{x, t) has the form 



^ l + expK-a^-3W2 3 
^ ' ^ 1 + exp(gi - 2ax) ' ^ ^ 

and we obtain the solution 

^ 1 - exp(-gi + 2ax + ln(2)) 

1 + exp(— gi + 2ax + ln(2)) + exp(mi — qi + ax — 3tips/2) 

For the new functions Q{x, t) — ui, M{x, t) — U2, we take the solutions of Eq. (2.34): 

'2exp(— gi) + exp(— 2aa;) + exp(mi — qi — ax — 3tips/2) 



H{x,t) = ln(^ 



U3 



1 + 2 exp(go — 2ax) 

+ exp(mo - qo- ax - 3tv?3/2)^ , 

/ , 2 + exp(gi)\ 
- 1 + exp Uo - ?i + 2ax + In — — — 

V 2 + exp(go)/ 

(. ( ^ , 2 + exp(fli) 

X 1 + exp Uo - ?i + 2ax + In — 

V V 2 + exp (go) 

( , , 1 exp(mi + go) + exp(mo + gi) ^^¥'3\VL.,x 
+ exp(^-g, + ax + ln ^^--^-^ 3— jj (2.41) 

By analyzing this formula, it is possible to understand the structure of the translation 
constant. Indeed, in the variables t\ — —ax + 3t<y93/2, T2 = ax + 3i</73/2, the solution 
(2.41) has the form 

-Cl exp(ri) + C2 exp(r2) 

^i3(Tl,T2) = 



1 + Cl cxp(ri) + C2 exp(r2) ' 



^ exp(gi)(2 + exp(go)) 
^ exp(go + mi) + exp(gi + mo) ' 
exp(go)(2 + exp(gi)) 

C2 — 



exp(go + mi) + exp(gi + mo) 



Lemma 2.1 The setu{x,t), M{x,t), Q{x,t) of solutions of the FhNS equation (2.34) 
of the form (2.38), (2.42) forms a semigroup with operation o (1-11) and the following 

properties: 

(a) commutativity, M{x, t) o Q{x, t) = Q{x, t) o M{x, t); 

(b) associativity, u{x, t) o {M{x, t) o Q{x, t)) = {u{x, t) o Q{x, t)) <> M{x, t); 

(c) any element is the unity for itself M{x, t) o M{x, t) — M{x, t). 

Proof Let us consider two solutions with different sets of constants Cj, Vi, i = 1,2 

n^Tiix^t). (2.43) 



Q{x,t) 
M{x,t) 



- exp(ri + Cl) + exp(r2 + C2) 
1 + exp(ri + Cl) + exp(r2 + C2) ' 

- exp(ri + Vi) + exp(r2 + V2) 
1 + exp(ri + Vi) + exp(T2 + V2) ' 
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By using Theorem 1.2, we find the solution 

u{x,t) = M{x,t)oQ{x,t) 

_ -exp(ri)(exp(Ci) + exp(Vi)) + exp(r2)(exp(C2) + exp{V2)) 
1 + exp(ri)(exp(Ci) + exp(yi)) + exp(r2)(exp(C2) + exp{V2)) ' 
T, = Ti[x,t\, (2.44) 

which imphes the properties given in Lemma 2.1. 

Note that the action of the semigroup is a nonhnear "time-translation" of the solution. 

As was noted in [2], p. 192, [3], p. 64, [9], and [12], it is known that Eq. (2.34) has 
solutions with singularities (monsters, contrast structures) whose role cannot be explained. 

We have the following theorem. 

Theorem 2.2 Suppose that the function Q{x,t) is a solution of the equation 

Qt - hQxx - ViQ + ^sQ^ = 0, (2.45) 
the function U{x, t) is a solution of the linear parabolic equation 

Ut + ^^ai/a^U, - y^s/a'f/.. = 0, (2.46) 
and the compatibility condition 

Q. + Qiai - i±)aQ - UJU) = (2.47) 
is satisfied. Then the function u{x,t) is a solution of Eq. ( p.34| ) and has the form 



Q{x, t)U{x, t) exp(t<^3(l + /(2a^))) 
^ Q{x,t)exp{aix) -U{x,t)exp{t<^3{l + ay{2a^)))' ^' ' 

Proof. To prove this theorem, we substitute the ansatz of the solution of the third 
mode 

, , exp{Hix,t))Mix,t) +exp(2Hix,t))Rix,t) -Qix.t) 

iilx t] = 

^ ' ^ 1 + Z{x,t)exp{H{x,t)) + exp{2H{x,t))Zi{x,t)+exp{3H{x,t))R{x,t)' 

(2.49) 

into Eq. (2.34), collect similar terms with equal powers of the exponentials, equate them 
with zero, and obtain a system of nine equations. As previously, from these equations 
we successively exclude the second derivatives of all functions contained there. Since the 
ansatz is invariant, there is a remarkable fact common for the algorithm of constructing 
the solution. Namely, just as in the theorems presented above, the equation becomes 
factored at some step. This equaton contains various hints, which allow us to correct the 
ansatz (2.50) and to pass to the next iteration. 
Equation (eg8) has the form 

2aV3 + ^a^Ht + ^fsHj - ^^H,^ = 0. (2.50) 

Its solution has the the form 

H{x, t) = - ln([/(a;, t)) + aix - {2a^ + al)t^-i / {2a^) , (2.51) 
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and thus implies (2.47). 

Equation (egO) is exactly Eq. (2.34) for the function Q{x, t). Hence the function Q{x, t) 
is a solution of this equation. 

Equation (egl) at the second iteration step implies the compatibility conditions. 

The second iteration of Eq. (eql) implies the compatibility condition in the form of the 
Riccaty equation, in which we study only the situation corresponding to the sign "plus" . 

Further, at each successive iteration, we have some versions dur to the fact that there 
situations not considered earlier. We choose 

Z{x,t) = -{M{x,t) + Q{x,tf)/Q{x,t), 
M{x,t) = ~Q{x,tf/2, 

R{x,t) = -{2Q{x,t)Zi{x,t) + Q{x,tf)/2. (2.52) 

By substituting all these expressions into the ansatz (2.50), we obtain the desired solution. 
The theorem is proved. 

There is the following analogy with the theory of Korteweg-de Vries equtaions. In this 
context, relations (2.48) can be considered as an analog of Miura transformation. 

Let us prove that all known solutions of the FhNS equation can be calculated in terms 
of the solutions of the a linear parabolic equation. We find the solution of the Riccary 
equation (2.48) and, choosing the upper sign, obtain 

Q(x, t) = ^(a:,t)exp(-aix) 

Co — /(a?7(x,t) exp(— aix) (ix 

The compatibility condtions impose strong restrictions. Hence, for example, we have the 
pair of functions 

Q(x,t) ^ l + ^iexp(2aa;) 



1 — Ci exp(2aa;) 

U{x,t) = Ci + exp(-2ax), ai = -a. (2.54) 
the substitution of which into (2.49) implies the following solution of Eq. (2.34): 

/ ,x 1 + Ciexp{2ax) 

Uix,t] = -; -. r -. -— . (2.55) 

— 1 + Ci exp(2ax) + exp(ax — 3tLp2,/2) 

And the solution of (2.47) 

U{x,t) = (l + 2aC2exp(—2ax)^ exp (j^a + al)x + {a — al)(a + al)t(p-i/{2a^)^, 

ai = IaV2. (2.56) 

gives solution of the equation (|2.34|) 

^ 2aC2 + exp{2ax) 

^'^ 2aC2-exp(2ax) + exp(aa;-3t(^3/2)' ^' ^ 

These solutions have a specific characteristic feature, which can either appear or disappear. 
Such solutions were considered in [3], p. 64, [9], [12]. In what follows, we consider 
Example 3 constructed by this method. In this example, two singularities simultaneously 
appear in the modified FhNS equation. 
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There is the following statement: to the WELL-KNOWN SOLUTION of Eq. (2.34) 

1 — exp(2ax) 

n = £_l i (2.58) 

1 + exp(2ax) + exp(ax — 31(^^/2) ' 

there corresponds the solution of the LINEAR PARABOLIC EQUATION (|]4|) 

U{x,t) = (— 1 + exp(2ax)) exp((— a + al)x + (a — al)(a + al)t(y93/(2a^)), 

ai = IaV2. (2.59) 

It turns out that the complex solutions of the linear parabolic equation can be recalculated 
into the solutions of the FhNS equation and vice versa. Therefore, the complex solutions 
of these equations must be studied. 

Let us consider the FitzHygh-Naguma-Semenov (FhNS) equation (2.34) with variable 
coefficients. 

By Theorem 2.1, for the first mode and separately for the third mode, there exist the 
following solutions. We have the following theorem. 
We have the following theorem. 

Theorem 2.3 Let = kQ{x,t), ipi = (pi{x,t), ip^ = (p^{x,t) be smooth functions. Sup- 
pose that the function M{x, t) is a solution of the equation 

Mt - koM,, - ipiM + ipsM^ = 0, (2.60) 

and the following compatibility condition imposed on its coefficients is satisfied: 

-V2 / ; — -— dx V^ + 6MV3V^-4Co(t)v"^ 

\J koy/rri2 / 

+QV2M,koip3 - V2Mko,ips + V2Mkoip3, = 0, (2.61) 
where 1712 = ko{x, t)ip^{x^ t) Then the function 

l-exp(/) J ^2koix,t) 

is a solution of the FhNS equation (2.34). 

Proof. The proof follows from Theorems 1.2 and 2.1 and can be carried out by a 
direct substitution of the ansatz of the first mode, taking into account the fact that the 
coefficients in the equation are variable and M = 0. One compatibility condition is 
identically satisfied. One compatibility condition remains. Hence the set of solutions is 
not empty. The following solution illustrates this theorem. 

For example, we set ko = 1, (pi = ip^. We also assume, that we want to dress 
the solution Q = 1. The compatibility conditon (2.62) turns into the equation for the 
function ipi. 



ifu - 3j2(piipi^ + {(pij /ipi + ifi^^ = 0. (2.63) 



We choose one of its solutions, for example, (fi = C2tanh.^{x/2J3C2V^)/S\/2 
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The solution of the FhNS equation 

Ut - u^x - - u^) = 0, (2.64) 

has the form 

') = /. = -p (^) cosh"' (^) . (2.65) 

This solution is real. If the constant is C2 > 0, then this solution describes wave interaction 
. At the initial time moment, the solution has a singularity. Then the solution becomes 
smooth. For negative values of this constant C2 < 0, this solution describes a periodic 
linear structure . 

Thus one can construct a solution of the equation with variable roots of the source-sink 
function. For example, if ipi is given to be a periodic function, then ip^ is constructed 
from solutions of Mathieu and Hill type equations. 



Theorem 2.4 Let ko = ko{x,t), (pi = !fi{x,t), ip^ = (fs{x,t) be smooth functions. Sup- 
pose that one of the compatibility conditions on the coefficients is satisfied, namely, the 
first compatibility equation 

QAx,t) = , (2.66) 

or the second compatibility condition 

^ I -2VkoUQ - 2UQ\, + 2^,QUx 

Q.(a:, t) = (2.67) 

Suppose that the function Q{x,t) is a solution of the FhNS equation 

Qt - koQx, - ^iQ + = 0, (2.68) 
and the function U{x, t) is a solution of the linear parabolic equation 

Ut + 2LpikoUx - koUxx ~ ^ ( / 'fitdx + Lpi + k^Lpl - k^Lpi^ = 0. (2.69) 
Then the function 

uix t) = Qix,t)U{x,t) 

^'^ Qix,t)exp{Jifidx)-Uix,ty ^ ' 



IS a so 



lution of the FhNS equation ( p.34|) . 



Proof. This statement can be proved by the substitution of the ansatz of the solution 
of the third mode ( p.49| ), just as in the proof of Theorem 2.2. This theorem is a partial 



case in the analysis of the one of the two branches, which are determined by the choice 
of the compatibility condition. 

The method has the following advantages: 

a) the algorithm is iterative; 
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b) at each iteration step we solve one ordinary differential equation of the first order 
or one algebraic equation (cf. [17], where one has to solve a system of ordinary differential 
equations); 

c) the method allows one to obtain new solutions; it is necessary only to know which 
mode has nontrivial solutions; to this end, there are some considerations the discussion 
of which is beyond the framework of this paper. 

Note that formulas (2.41), (2.42) can be also obtained in a different way, while formulas 
(2.49), (2.63), (2.72) are constructed by this method. 

EXAMPLE 2. Let us consider a modified FitzHygh-Nagumo-Semenov equation, 
which is a special case of the equation ( |1.1|) . In this case there also exist a lotof so- 
lutions, but we write only one solution correspondign to the second mode, since it differs 
from all other solutions. The change ( p.. 21) reduces our equation to an equation the fourth 
order of homogeneity. 



Suppose that Eq. (|TTl|) has the form 

{ai — hiu)ut — {ai — hiu)uxx — hi{ai+hiu)ux—2hi{uxY—ai{a\+h\)u{l—u^)/2 = 0, (2.71) 
The solution has the form 

1 — expfaix) 

y r f) = 2 72) 

l+exp{aix)+exp{{{-hj-3al)t + 2{hi + ai)x)/A) ^ ' 

and describes the development of the second wave from small perturbations. 



This solution can be used for modeling phenomena related to phase transitions [10 



One of the elementary equations for which this effect is preserved for ai = has the form 

Ut-UxX + hiUx + 2{uxy/u = 0, (2.73) 



EXAMPLE 3. We consider the modified FhNS equation ([Ol) : 

ut + 2{hifuxJ^-i - Hi{l + Au)ux - ^3u{l - u^) = 0, (2.74) 
The solution has the form 

^'^ l-exp(x(^3/(2/ii))+exp((M + a;)(^3Mi) 

This exact solution shows how a solution of the problem with a smooth initial condition 
of double singularity is developed during a finite time interval. A discontinuity arises at 
t = 1.4. This discontinuity is of a special structure such that the left-hand discontinuity 
does not move, while the right-hand discontinuity moves to the right. 

Note that if the sign of the parameter ip^ is changed (i.e., we have cp^ = 1), then 
Eq. ( |1.1[ ) becomes an inverse parabolic equation. However, the solution ( p..lD exists and 
is bounded. Small perturbations of the initial condition relax in time. 

EXAMPLE 4- We consider the solution correspondign to the second mode for the 
equation the fourth order of homogeneity. Suppose that the modified FhNS equation 
(|1.1|) has the form 

(1 - u)ut - (1 - u)uxx + (-l/oi + u/ai- 2aiU(f3)ux - 2{ux f - (f3u{l - u^) = 0. (2.76) 
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This equation has a solution of the form 



t] = , 1-^^P(^/^1) ^. (2.77) 

1 + exp(x/al) + exp(x/al — tip^) 

The solution shows how the structure develops from an initial state whose range is 
[—0.5, 1] to a traveling wave u G [0, 1]. Note that it is possible to construct many solutions 
of such type, i.e., solutions that describe how a certain wave evolves from some intial state. 

EXAMPLE 5. We consider the modified Fisher-Kolmogorov-Petrovskii-Peskunov 
equation 

ut - koUj^j, - {-2kof2/h2 + h2u)uncx - y^iu + ip2u'^ = 0. (2.78) 

We consider an important example of solutions to the Fisher-Kolmogorov-Petrovskii- 
Piskunov modified equation and show that its solutions can be expressed in terms of higher 
transcendental functions such as Bessel functions or hypergeometric functions. We also 
show how they are related to the Backlund transformation for the Biirgers equation. 

Theorem 2.5 Suppose bi = 0, (f^ = 0, 

hi = -2koifi2/h2 (2.79) 
in equation (pT3|). The function U{x,t) is a solution of the linear parabolic equation 

Ut = r{x,t)U + koU,,, (2.80) 

with the potential 

r(x, t) = exp(2x//i2 + Ci(t)). (2.81) 
The function U{x, t) has the form 

U{x, t) = v{z{x, t)), Z{x, t) = exp (x/{2h2) + Ci(t)/2)) (2.82) 

and the function v{z) is a solution of the equation 

2hlzip2v{z) + hjifiiv'iz) + 2ko(p2v'{z) + 2kQZLp2v"{z) (2.83) 

Then the function u{x,t) is a solution of Eq. ( ^.78] ) and has the form 

u{x, t) = {2ko/h2)dMUix, t))), (2.84) 

Proof. We substitute ( p.84|) into the original equation (|2.78| ) and obtain a nonlinear 
homogeneous equation from which we exclude the derivatives Ut, U^t with the help of 
Eq. ( p.80|) . Then, for the parameters specified in the assumptions of the theorem by 
the change ( 2.83 ), we obtain ordinary defferential equations for higher transcendental 
functions. The proof of Theorem 2.5 is complete. □ 

It should be noted that the change ( |2.84D coincides in form with the Cole-Hopf change 
of variables for the Biirgers equation. 

The author is grateful to V. G. Danilov and S. Yu. Dobrohotov for constant attention 
to his work and useful discussions and to V. P. Maslov, and A. D. Polynin for constructive 
advice. 
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